Simulations of energetic beam deposition: from picoseconds to seconds 



Joachim Jacobsen^, B.H. Cooper, and James P. Scthna 
Laboratory of Atomic and Solid State Physics and Cornell Center for Materials Research, Cornell University, Ithaca, NY 

14853-2501, USA 

^Present address: Haldor Topsoe A/S, Nymollevej 55, DK-2800 Lynghy, Denmark. 



(February 1, 2008) 



X 
S3 



We present a new method for simulating crystal growth by 
energetic beam deposition. The method combines a Kinetic 
Monte-Carlo simulation for the thermal surface diffusion with 
a small scale molecular dynamics simulation of every single de- 
position event. We have implemented the method using the 
effective medium theory as a model potential for the atomic 
interactions, and present simulations for Ag/Ag(lll) and 
Pt/Pt(lll) for incoming energies up to 35 eV. The method 
is capable of following the growth of several monolayers at 
realistic growth rates of 1 monolayer per second, correctly ac- 
counting for both energy-induced atomic mobility and ther- 
mal surface diffusion. We find that the energy influences is- 
land and step densities and can induce layer-by-layer growth. 
We find an optimal energy for layer-by-layer growth (25 eV 
for Ag), which correlates with where the net impact-induced 
downward interlayer transport is at a maximum. A high step 
density is needed for energy induced layer-by-layer growth, 
hence the effect dies away at increased temperatures, where 
thermal surface diffusion reduces the step density. As part of 
the development of the method, we present molecular dynam- 
ics simulations of single atom-surface collisions on flat parts 
of the surface and near straight steps, we identify microscopic 
mechanisms by which the energy influences the growth, and 
we discuss the nature of the energy-induced atomic mobility. 

PACS numbers: 81.15.Jj, 68.55.-a, 61.72.Cc 



I. INTRODUCTION 

The overall goal of being able to manufacture nanoscale 
devices while controlling chemical composition and crys- 
tal structure has been the driving force behind a tremen- 
dous effort in numerous research groups over the past 
decades, and the ability to grow crystal surfaces in a 
layer-by-layer fashion has been a key issue The use 
of energetic particles has offered a promising possibility 
of gaining control of the growth process . 

There are numerous examples where energetic beams 
have been used to control and improve the properties 
of the grown materials. Ion Beam Assisted Deposition 
Q, Ion Beam Sputter Deposition js) and Ion Beam Di- 
rect Deposition have been used to lower the epi- 
taxial growth temperature and improve smoothness. In- 
creased control of interfacial roughness has been achieved 
through sputter deposition 0], and Giant Magnetoresis- 



tance has been improved both using both sputter depo- 
sition [p|,p^ and using direct ion beam deposition [pd]. 

Not surprisingly, there can be an optimal window for 
the energy per incoming particle: a certain energy is 
needed to increase atomic mobility at the surface, but 
too high an energy can cause a drastic increase in defect 
formation. There is a growing literature on experiments 
which map out parameter space to find optimal values 
for various systems, and the optimal energy window is 
often found at relatively low energies. For example, Ra- 
balais et al. found that for silicon ion-beam epitaxy very 
smooth growth was obtained using 20 eV particles 0. 

There is an obvious need for improving our un- 
derstanding of how the energy influences the growth. 
What are the microscopic mechanisms by which the en- 
ergy changes the growth morphology? Recent experi- 
ments have illustrated how the energetic ion-surface col- 
lisions can influence submonolayer island densities during 



growth of a Pt(lll) surface by sputter deposition [n2 13 



or Ion Beam assisted deposition ]lj] . Several experimen- 
tal groups are operating or building equipment for 
Ion Beam Direct Deposition with high control of beam 
energy and angle of incidence [ p^ , and we can expect 
much more detailed experimental information on the ef- 
fect of the energy in the near future. 

Molecular dynamics (MD) simulations using model po- 
tentials, from which very detailed information on micro- 
scopic mechanisms can be obtained, have proven to be a 
very useful tool in studying energetic ion-surface ll^-|2l|] 
and cluster-surface collisions. There are several 
questions that can be directly addressed using MD. What 
kind of atomic rearrangements occur in the ion-surface 
collisions as a function of energy? Is ballistic motion 
or local heating the right picture of the energy-induced 
mobility? What are the time and length scales of the in- 
duced mobility? For example, Villarba and Jonsson |l^] 
studied low energy (10 and 20 eV) impact of Pt atoms on 
a Pt(lll) surface, and identified push-out events where 
atoms impinging close to descending steps are incorpo- 
rated into the growing layer. The net effect is that the 
step edge grows horizontally favoring a layer-by-layer 
growth mode. They found that the non-thermal effect 
of the collision was over in a few picoseconds, and the 
ranges of the collision-induced atomic rearrangements 
were a few lattice sites. 

One issue, which MD fails to address, and which is cru- 
cial in understanding most experiments using energetic 



1 



particles, is what is the relative importance of collision- 
induced atomic mobility on the one hand, and the ther- 
mal surface diffusion on the other? MD simulations at- 
tempting to simulate the entire growth process must use 
deposition rates on the order of 10^*^ monolayers per sec 
(ML/s) |l^[20[|, and neglect the effect thermal surface 
diffusion on longer time-scales. Furthermore, subsequent 
collision events might affect each other directly in ways 
they would not at the much lower experimental deposi- 
tion rates (typically 0.01 - 1 ML/s). If we want to ad- 
dress this question of the relative importance of collision- 
induced mobility and the thermal surface diffusion, we 
are faced with a tremendous problem of time-scales. 
Each atom-surface collision is a pico-second event, to re- 
solve which you need femto-second numerical time steps. 
At the other extreme, in the typical experimental situ- 
ation there is one such collision per lattice site per 1 to 
100 seconds. 

In this paper we present results of a new simulation 
method capable of overcoming this 15 orders of magni- 
tude gap in time-scales. We have combined two well 
known simulation techniques in a new way. We do a tra- 
ditional Kinetic Monte-Carlo (KMC) simulation of the 
thermal surface diffusion in between the rare collisions, 
and for each of these collisions we do a small length-scale, 
short time-scale molecular dynamics simulation. We call 
this a hybrid Kinetic Monte-Carlo Molecular Dynamics 
(KMC-MD) simulation. We will show the results of de- 
positing several monolayers at IML/sec, while treating 
every collision explicitly in MD simulations. As model 
systems, we have chosen Ag(lll) and Pt(lll) homoepi- 
taxy, and we use Effective Medium Theory (EMT) |2|] 
as a model potential. An idea similar to KMC-MD has 
been used to simulate damage production in ion implan- 
tation of silicon . 

We use effective medium theory, which is known to 
give a good qualitative and to some extent quantitative 
description of these metals [ p3p^ ]. We do not expect 
these simulations to accurately predict all features of 
the growth by energetic beam deposition of Ag(lll) and 
Pt(lll) However, we do hope to gain useful in- 

sight into how the use of energetic particles can influence 
growth, microscopic mechanisms, energy-induced atomic 
mobility, and the interplay between the energy-induced 
mobility and thermal surface diffusion. 

The paper is organized as follows: In section |l^, we de- 
scribe the KMC-MD method in detail, going through 
how the KMC part of the simulation is done (sec- 
tion II A), how the MD part is done (section II C), how 



to go from the d iscre te KMC to the conti nuou s MD sim- 
ulation (section IIB), and back (section II D). The re- 
sults section 



part 



III 



is divided into two main parts 
- MD simulati ons of single energetic atom-surface col- 
lisions (section III A ) and KM C-MD simulations o f the 
entire growth process (section IIIB). In section III A 



avoid unphysical reflections of the supersonic shock-wave 
induced by the impact. We then go through some im- 
portant energy induced microscopic mechanisms we find 
for Ag— >Ag(lll) and Pt— +Pt(lll) energetic impacts, fo- 
cusing on the energy-induced upward and downward in- 
terlayer mobility at straight step edges. We round up 
giving a short s ummary and discussion of these mecha- 
nisms in section III A q. In section III B , we present our 



we start our discussion of the choice of MD system size 
(section [II A 1 ) , where we show that large MD systems 
and Langevin damping on the boundaries are essential to 



KMC-MD results for the submonolayer structure and for 
the surface roughness after the deposition of a few atomic 
layers; first for Ag/Ag(lll) and then for Pt/Pt(lll). We 
find the smoothest growth when both the net downward 
interlayer mobility and the step edge density are large. 
In section we give a concluding discussion of our re- 
sults for simulations of growth by energetic beam depo- 
sition, and in section ^ we give a general discussion of 
the KMC-MD method. 



II. THE KMC-MD METHOD 

In the introduction above, we stated that a main dif- 
ficulty in simulating crystal growth by energetic deposi- 
tions is the tremendous gap in time scales between the 
picosecond collisions and the deposition, typically slower 
than one impact per lattice site per second. The idea 
behind the KMC-MD method is precisely that the de- 
positions are rare events. For most of the time during 
the crystal growth, there is no non-thermal atomic mo- 
bility on the surface. Only the thermal diffusion of atoms 
is active. Kinetic Monte-Carlo is a very efficient way of 
evolving the surface for this time in between the energetic 
collisions, given a model of the thermal diffusion. 

Once in a great while (on the time scale of diffusion) 
an energetic collision occurs. These collisions may dis- 
place the incoming atom from the impact site, as well 
as rearrange the surface atoms. But for homoepitaxial 
metal systems, the energy transfer is very efficient, and 
the incoming atom generally slows down very quickly. 
Also, in the crystalline environment, the excess energy 
dissipates away very quickly due to phonons. Conse- 
quently, non-thermal mobility is limited to a very short 
time, te, after the collision. Furthermore, for not too 
high incoming energies, the collision-induced mobility is 
fairly short ranged. This makes MD of the collisions a 
suitable method. 

In the KMC-MD method, we evolve time using a KMC 
simulation in a standard way |2^,|2^,^,^. Diffusion 
processes happen sequentially according to their relative 
rates and, with a probability proportional to the depo- 
sition rate, new atoms are introduced above the surface 
with a specified kinetic energy and angle and aimed at 
a random point on the surface. When this happens, we 
set up an MD simulation to correctly describe the atom- 
surface collision. We include only the local region in the 
vicinity of the impact site in the MD simulation, and run 
for only a short time. We then feed the end result back 
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into the KMC simulation, and continue. In sections y 
A-D below, we describe the details of the different parts 
of the KMC-MD method. 

We model the atomic interactions using effective 
medium theory — in the MD simulations of the 
energetic collisions we use the EMT forces, and in the 
KMC~simulations we use a comprehensive set of EMT 
energy barriers for various atomic diffusion processes. 

A. The KMC lattice simulation 

Kinetic Monte-Carlo has become a 
standard method for doing lattice simulations of crys- 
tal growth [^^5|j29|j2^ . In short, we make a complete 
table of the active atomic diffusion processes at the given 
temperature. For each type of diffusion process we evalu- 
ate its rate as = z/exp (—Ei/ksT), where we assume a 
common prefactor oi v — lO^^s"^ for simplicity. Table || 
lists the processes we include together with their EMT 
energy barriers Ei for the two systems (Ag/Ag(lll) and 
Pt/Pt(lll)). The surface atomic configuration is spec- 
ified by the occupancy by atoms of the fcc(lll) lattice 
sites. For every atom on the surface we examine if it can 
potentially make a lattice jump in any of the processes 
from table |, and if so we add these specific atomic jumps, 
and their rate to a list of all the possible potential diffu- 
sion processes the given surface configuration can evolve 
by. Also included in the list is the deposition of a new 
atom — with a rate given by the deposition flux times the 
surface area. Now, in every loop of the program, with a 
probability proportional to its rate, one particular atomic 
jump is chosen from the list of potential processes. The 
surface configuration and the list of potential processes 
is then updated accordingly, taking advantage of the fact 
that there will be only local changes. Because there is 
a lattice jump or a new deposition in every loop of the 
program, KMC is a very efficient method, and can bridge 
one gap of time-scales: the gap between fast thermal dif- 
fusion processes (10^ to 10^ per sec.) and the deposition 
rate. We depart from the standard way of doing crys- 
tal growth with KMC in the way we handle deposition 
events. Usually in KMC simulations when a deposition 
event is chosen to happen, a new atom is introduced at a 
random surface lattice site. Instead we do something new 
— a complete MD simulation of the atom-surface colli- 
sion with the specified energy and angle and an impact 
parameter for the incoming atom chosen at random. 

In simulations of thermal crystal growth and relaxation 
the KMC method has been used in two ways: the energy 
barriers giving the rates can be obtained from a model 
potential energy, or they can be used as adjustable pa- 
rameters to fit growth experiments or to study the effect 
of a particular process on the resulting surface morphol- 
ogy. In our case, to make the two parts of the simulation 
consistent with one another, the KMC energy barriers 
must be obtained from the potential energy used in the 



MD simulations. 



B. From discrete KMC to continuous MD 

When the KMC algorithm chooses a deposition event, 
a random point in a plane above the surface is picked. 
From this point, a straight line at the ion beam angle 
is then followed until it crosses a horizontal (lll)-plane 
near (within a surface unit cell of) a lattice site occupied 
by an atom, which we will call the impact atom. With 
respect to this expected impact atom, the cluster is set up 
based on the occupancy of sites in the KMC-simulation, 
as discussed below. 
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TABLE I. Effective Medium energy barriers in meV used 
in the KMC simulations. In detachment processes an atom 
stays at the surface, but dissociates from in-layer neighbors. 
Edge diffusion is for atoms moving along island edges. Ni 
and Nf are the initial and final in-layer coordinations of the 
moving atom. The transition state is labeled A or B ii the 
motion is along a (100) or (111) micro-facet. Dimers can 
diffuse via a single atom mechanism in which one atom moves 
along the edge of the other. The atomic moves are illustrated 
in detail in refs. [Q. The " — " for Ag indicates that the 
barrier is not distinguished from the one immediately above. 
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Before introducing the energetic atom, the system is 
equihbrated at the specified temperature. We find that 
an efficient way to do this is to give every atom a veloc- 
ity picked from a Maxwell-Boltzmann distribution cor- 
responding to twice the temperature. Since the atoms 
start out on perfect lattice sites, this gets the total en- 
ergy per atom roughly right. We then need only to run 
for 0.27 ps to get a reasonably equilibrated system. Fi- 
nally, the atom is placed three horizontal (fff) planes 
upward along the straight line introduced above, and is 
started with the specified kinetic energy. 



C. The MD continuous simulations 
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FIG. 1. Cluster setup used in the MD simulations. The 
MD-atoms (white) have zero Langevin friction coefficients, 
the Langevin atoms (gray) have positive Langevin friction 
coefficients, and the static atoms (black) are immobile and 
fixed at lattice positions. Left: view from glancing angle, 
and view of a vertical cut. Right: view from below of the 
3 types of atoms. We include atoms in shells bounded by 
(111) and (100) planes at specified distances to the impact 
atom (marked with x). Two integers specify the region of 
MD-atoms atoms: Ni is the depth in number of horizontal 
(111) planes and N2 is lateral size in number of (111) or (100) 
planes, as shown. The region of Langevin atoms is the sur- 
rounding N3 number of (111) or (100) planes, and likewise 
the region of static atoms is the surrounding N4 number of 
(111) or (100) planes. The shown setup has Ni = 5, N2 = 6, 
N3 = 3, and N4 = 4, or N = (5, 6, 3, 4). (Ni = N2 = 1 would 
give exactly 1 MD-atom). 



Following DePristo and Metiu we set up a clus- 
ter consisting of three types of atoms. The atoms in the 
immediate vicinity of the impact site evolve classically 
according to the EMT forces only (we call these MD- 
atoms). Surrounding the MD-atoms we have a shell of 
Langevin atoms, which are subject to the EMT forces, 
and in addition to a friction force and a randomly fluc- 
tuating force. To fix the geometry a shell of static atoms 



surrounds the shell of Langevin atoms (Figure |lj). The 
equations of motion for the dynamic atoms are: 

= -Vr^VEMT{{r}) - m,6^ + h{t) (2.1) 

where and mi are the three Cartesian coordinates and 
the mass of atom i, Vemt{{t}) is the EMT potential 
energy as a function of all atomic coordinates {r}, is 
the Langevin friction coefficient of the atom i (the use of 
the index will be modified below), t is time and f{t) is 
the fluctuating force obeying the fluctuation-dissipation 
theorem: 



{mum = 2^,m,kT5{t). 



(2.2) 



The central MD-atoms h ave zero friction coefficient, = 
0. We integrate equation 2T in time using the algorithm 
proposed by Allen and Tildesley |Q, which reduces to 
the Verlet algorithm for = 0. 

Having a shell of Langevin atoms which surrounds the 
MD-atoms serves several purposes. It allows us to equili- 
brate the system at a specified temperature before shoot- 
ing in the energetic atom; it mimics the contact of the 
MD-atoms with an infinite heat bath (the crystal), ensur- 
ing that the deposited energy doesn't permanently heat 
up the system; and furthermore, as we shall see later, it 
allows us to use a smaller system size without unphys- 
ical refiections of energy from the boundaries affecting 
the atomic motion near the impact site. We will re turn 
to a discussion of this latter point in section |III A . We 
always use — Z and A'4 = 4, except for a test cal- 
culation with no Langevin atoms {Nj, — 0). Having 3 
atomic planes of Langevin atoms, we use the possibility 
of having different Langevin coefficients for each plane. 
We label these coefficients ^ = (■^17^2,^3), where ^1 is for 
the plane neighboring the MD-atoms, and ^3 is for the 
plane neighboring the static atoms (fig. |l|). 

Apart from how the cluster boundaries affect the en- 
ergy dissipation away from the impact site, another cri- 
terion for choosing the size of a cluster simulation is the 
range of the transient mobility induced by the collisions. 
In the energetic collisions, atoms are usually displaced a 
few lattice sites from the point of impact, and it is im- 
portant that these atoms stay in contact with the region 
of MD-atoms. Occasionally, incoming atoms do not stick 
to the surface, but are reflected back. For rare impact 
parameters on step edges, incoming atoms may travel 
rather long distances along the surface before sticking 
and thermalizing. In long simulations with thousands of 
depositions such events will be encountered. However, at 
the end of an MD simulation where the lattice sites of 
the atoms have to be identified, we check if any dynamic 
atoms left the physical region of the cluster calculation. 
Our KMC-MD system sizes were sufficiently large that 
this happened in fewer than 2% of the deposition events. 

Every MD simulation is followed for 5 ps. After this 
time, the average kinetic energy per atom in the local 
impact region is once again ~ 3/2 /cT, and all further 
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atomic mobility will be thermally activated and well- 
described by the KMC lattice simulation. 

We find it important to use a variable time step in 
the MD simulations of the collisions. However, based 
on Figures || and we choose to do this in the simplest 
possible way. The very high kinetic energies are all found 
within the first 0.5 ps of the collisions, and after this time 
the atoms have much more moderate velocities. We use 
the Verlet algorithm, with a time step of dti = 1.08 fs 
for the first 0.5 ps, and subsequently we use dt2 = 5.4 fs 
for the remaining 4.5 ps. By comparing to much shorter 
time steps, we have checked that this scheme, for energies 
up to 30 eV, is sufficient to correctly identify the lattice 
sites of all atoms after 5 ps, whereas using dt2 for the 
entire collision would give inaccurate results. 



D. From continuous MD to discrete KMC 

After having the done the MD simulation for the lo- 
cal vicinity of the impact, when all atomic mobility is 
once again thermal, we need to map the end result of 
the continuous simulation back to the lattice description 
of the KMC part of the simulation. In the thermal sur- 
face configuration every atom is oscillating around a local 
minimum of the EMT potential energy surface. To iden- 
tify the binding site of each atom, we simply minimize 
the energy with respect to the atomic coordinates. This 
uniquely identifies an fee lattice site for almost all the 
atoms. 

However, there is one exception to this. The presence 
of hep binding sites on the fcc(lll) surface give rise to 
off-lattice local energy minima which can trap atoms. In 
EMT, an isolated adatom at the surface has essentially 
the same energy when at an fee and an hep site, and has 
roughly equal probabilities of occupying the two types of 
sites. As a result, during growth islands can nucleate on 
hep sites as well as fee sites, which would possibly make 
the coalescence of each new layer a very complicated pro- 
cess involving surface dislocations separating regions of 
hep and fee stacking sequences. However, low tempera- 
ture growth experiments for Ag/Ag(lll) and Pt/Pt(lll) 
do not show the complicated behavior expected if islands 
of substantial size had significant probabilities of being 
on hep sites — for Ag/Ag(lll) the growth is known to 
proceed on the fee sites . More accurate total energy 
calculations for Pt/Pt(lll) |Q finds a signifieantly in- 
creased energy for the adatom at the hep site, suggesting 
that the strong binding at hep sites is an artifact of the 
effective medium theory, due to the lack of dependence 
of the energy on relative bond angles. For these reasons, 
and because it would make the KMC simulation infeasi- 
ble, we do not allow atoms to occupy hep sites. 

Since the hep sites are unfavorable, and thermal dif- 
fusion will rapidly shift atoms from hep metastable sites 
to neighboring fee sites, we need a procedure which takes 
the atoms on hep sites in the energy minimized config- 



uration, and puts them on one of the three neighboring 
fee sites. When doing this the important thing is not to 
break any bonds if, for example, a dimer on hep sites 
has to be displaced. For each atom that ended up on an 
hep site after the energy minimization, we move it to the 
neighboring vacant fee site of highest in-layer coordina- 
tion, or a random choice among these sites if there are 
several equivalent ones. 

With all atoms occupying fee sites, we can continue 
the KMC simulation of the thermal surface diffusion. 



III. RESULTS 

We present the results of our simulations in two sub- 
sections. Before showing the results for the growt h sim- 
ulations with the KMC-MD method in section IIIB , 
we present and discuss M D sim ulations of single atom- 
surface collisions in section III A below. We start out dis- 



cussing the choice of system size in the MD simulations, 
and how for small systems this affects the outcomes of 
the collisions. We then move on to classify and quan- 
tify the collision-induced processes influencing atomic 
mobility found at relatively low energies (< SbeV) for 
Ag— >Ag(lll) and Pt^Pt(lll). Unless otherwise spec- 
ified, the collisions occur at normal incidence and with 
random impact parameters. 



A. MD of single atom— surface collisions 

1. Choice of MD system size: avoiding reflected energy from 
the boundaries 

In atom-surface collisions, the deposited energy is dis- 
sipated away from the local impact region by shock- 
waves and phonons. Clearly, how much energy stays 
around for how long a time is crucial for how much 
atomic mobility the collision causes. When doing a clus- 
ter simulation of the collision, it is possible that the finite 
system size will affect the energy density in the impact 
region after the collision — outgoing energy waves will 
be reflected at the boundaries and come back to the lo- 
cal impact region. This is illustrated in figure ^, which 
shows results for 25 eV Ag— >Ag(lll) impacts for a sys- 
tem setup given by iV = (3,4,3,4) and Langevin co- 
efficients C = (0.005,0.010,0.015). The lower panel in 
the figure shows the kinetic energy averaged over many 
impacts for individual atoms near the impact site as a 
function of time. For the first 0.4 ps, this graph is in- 
distinguishable from the corresponding graph for a much 
bigger system. At 0.4 ps, we see that the second-layer 
atoms (#3) suddenly accelerate, and immediately after 
that the impact atom in the first layer (#2) gets a boost 
of kinetic energy. This sudden acceleration after 0.4 ps is 
absent in a calculation for a sufficiently big system, and is 
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due to energy reflected by the static atoms at the bound- 
aries of the cluster. The upper panel in the same figure 
shows the height of the impact atom Z2 above the first 
surface-layer as a function of time, for thirteen individual 
impacts. We see that this atom sometimes stays in the 
same layer, sometimes pops out of the surface layer and 
equilibrates as an adatom, and once ends up one atomic 
layer down. The fact that the impact atom can pop out is 
an energy-induced effect. It does not happen for thermal 
impacts. However, the figure illustrates that for system 
sizes this small, the impact atom gets a sizeable boost of 
energy reflected from the boundaries right at the criti- 
cal time: just where some trajectories take the atom out 
of the surface layer and some don't. And, as we shall 
see shortly, the probability of surface atoms popping out 
in the energetic impacts depends on the system size for 
small clusters. 
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FIG. 2. False energy reflections from the MD simu- 
lation boundaries: 25 eV Ag^Ag(lll) impacts. The lower 
panel shows the average over many impacts of the kinetic en- 
ergy of individual atoms at various positions as a function of 
time after the impact. The inset schematically shows the dif- 
ferent types of positions. 1 is the incoming atom and 2 is the 
impact atom. In the fee lattice, there are really 3 atoms in 
the position labeled 3, and 6 atoms in the position labeled 4. 
The system has iV = (3, 4, 3, 4) and ^ = (0.005, 0.010, 0.015). 
(Same as Si in Fig.^). The reflected energy from the bound- 
aries is evident. After an initial peak around 0.1 ps the atoms 
slow down. Then around 0.4 ps, atoms at position 3 (dot- 
ted curves) speed up, and immediately after that the impact 
atom accelerates wildly. The upper panel shows the vertical 
position Z2 of the impact atom for thirteen individual im- 
pacts. Starting out in the surface layer at Z2 = 0, the impact 
atom gets kicked towards the surface, where it bounces off the 
second layer and moves back out. Around 0.5 ps its fate is 
determined. It either stays in the surface layer, or pops out 
by one atomic plane (z2 — > 4); in a single case it squeezes into 
the second layer (z2 —4). The reflected energy assists the 
impact atom in popping out. 
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FIG. 3. System size and Langevin reduction of en- 
ergy reflection artifact. Average kinetic energy of the 11 
local atoms at positions 1 to 4 in Figure |^ for one hundred 
25 eV Ag — » Ag(lll) impacts in 4 different system setups. 
The initial surface temperature is 17 meV (200 K). The inset 
(same units as the main plot) is the same plot on a different 
scale, showing the initial drop in kinetic energy as the 25 eV is 
distributed over many degrees of freedom. The main plot fo- 
cuses on a bump in the local kinetic energy observed after 0.5 
ps, and shows the system size dependence of this bump. For 
small systems, the bump is too high due to energy reflected 
from the boundaries. The setup S2 has no Langevin atoms. 
The setup S3 has the same size as S2, but S3 has 3 shells of 
Langevin atoms, and the damping clearly reduces the amount 
of reflected energy. 

Figure || shows the kinetic energy of the 11 atoms 
neighboring the impact site (average of the 11 curves in 
the lower panel of figure ||), again averaged over many 
impacts, and now for four different system sizes and 
Langevin coefficients. This plot shows a bump in the 
local kinetic energy around 0.5 ps after the impact, and 
that the size of this bump depends on the system size 
and on the Langevin coefficients of the boundary atoms. 
The small system setup 5*1 is the same as in figure ^ 
and has a total of 757 dynamic atoms (MD-atoms -|- 
Langevin atoms), S2 and 6*3 both have 1093 dynamic 
atoms (for 5*3 some of these are Langevin atoms), and 
S'4 has 2741 dynamic atoms. As noted above, curves for 
smaller systems can deviate from those for bigger system 
sizes because of refiections from the boundary. That the 
curve for 5*1 deviates significantly from S'4 after 0.36 ps, 
whereas S2 and 6*3 deviate from S'4 only after 0.45 ps, 
reflects the fact that that the bigger the system the later 
the traveling reflected energy wave will affect the local 
impact region. We believe that the bump in the curve 
for S4 will not disappear as the system size is increased 
further; extrapolating from the times at which the curves 
for Si , S2 , and S3 deviate from S4 , we estimate that the 
reflected energy for S4 does not return until after roughly 
0.5 ps, at which time this bump has already reached its 
peak. The size of the bump is reduced from Si to S2 
by increasing the system size. It is then reduced fur- 
ther from S2 to S3 merely by tuning the Langevin coef- 
ficients on the boundary atoms. We have done an exten- 
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sive numerical exploration to find the optimal choice of 
Langevin coefficients, and found = 0.005, ^2 = 0.010, 
and ^3 — 0.15 to be a good choice. A Langevin coef- 
ficient too high and too close to the impact region will 
cause energy reflection at early times directly from these 
Langevin atoms. The value ^3 = 0.15 seems to be the op- 
timum for the outer most layer of Langevin atoms. We 
have also tried random Langevin coefficients within each 
layer, to avoid focusing of the reflected energy, but with- 
out finding improvements. It is evident that the curve for 
^3 is not completely converged: the local kinetic energy 
is affected by the boundaries for the system S3. However, 
we shall argue below that it is sufficiently converged that 
it does not significantly affect the atomic rearrangements 
in which we are interested. 

We note that the time at which the reflected energy 
comes back tells us that the energy wave travels at a 
supersonic speed. This is due to anharmonicities in 
the potential energy at the relatively large atomic dis- 
placements: the atomic collisions at the earliest times 
are presumably exploring the hard core of the poten- 
tial and hence can propagate faster than harmonic sound 
waves. Potential energy anharmonicities (e.g., backscat- 
tering from hard collisions) can also cause reflection of 
outgoing energy even from the MD-atoms following the 
true dynamics. We believe these anharmonicities cause 
the bump in local kinetic energy found around 0.5 ps in 
figure ||for the largest system (curve 54). 



2. Adatom/vacancy formation 




FIG. 4. Possible outcomes of 35 eV impact of Ag 
Ag(lll) initially at 200 K. The top row shows initial config- 
urations with atom 1 (the incoming atom) above the surface, 
and the bottom row the corresponding final configuration af- 
ter thermalization. Atoms shifting lattice site are numbered. 
The three right columns show examples of adatom/vacancy 
formation, while the left column is an example of an exchange 
event with no net effect on the substrate. The processes shown 
in the first and second column each account for about 30% of 
the impacts at 35 eV. The process in the third column is an 
exchange event in which four neighbors of the atom marked 
"x" forming a chain in the first and second layers each shift 
positions by one atomic distance. (Unlike earlier figures, here 
"x" does not mark the impact atom, which in fact occupies 
the site where the vacancy is found after the collision.) At 
25 eV, the process in the first column is the single most com- 
mon event; and the process shown in the second column is 
the most frequent way of forming an adatom/vacancy pair. 
Projected trajectories for this process are shown in Figure H. 



Figure ^ shows possible atomic rearrangements that 
result from 35 eV Ag Ag(lll) impacts on an initially 
flat surface. In one case (first column from left), there 
is an exchange process, in which the incoming atom gets 
incorporated into the surface layer, and a surface atom 
pops out. The impact has no net effect on the substrate. 
The only net effect of the incoming energy is that the re- 
sulting adatom ends up a couple of atomic spacings away 
from the impact site. It is an example of energy-induced 
short ranged horizontal mobility, which is very common 
for these low energy impacts. In the second column from 
the left, something more happens. The incoming atom 
gets incorporated into the surface, but two surface atoms 
pop out. Compared to a thermal deposition, which would 
result in one adatom on the surface, the energetic impact 
causes the formation of an additional adatom/vacancy 
pair. This can potentially have a significant infiuence on 
the subsequent growth. Since these two adatoms are cre- 
ated very close to each other, they have a large likelihood 
of meeting and f ormin g a dimer. We shall return to this 
point in section III B| . Figure ^ also shows an example 
in which two surface vacancies are formed, and a total 
of three atoms end up outside the surface layer. Two of 
these atoms have met to form a dimer during the collision 
process. 
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Before we turn to a discussion of the probability of 
making these adatom/ vacancy pairs as a function of the 
incoming energy, let us look into the mechanism by which 
they are formed. Figure ^ shows the atomic trajectories 
for the impact shown in figure ^ column 2, in a vertical 
cut. This process is the most frequent way of forming 
adatom/ vacancy pairs at 25 and 35 eV. We see atom 
1 coming in, squeezing atoms 2 and 4 down and apart. 
After bouncing off the atoms in the second surface layer, 
the atoms 2 and 4 push atoms 3 and 5 out of the surface. 
Atom 1 takes the place initially occupied by atom 4. The 
picture that emerges is very ballistic. A collision sequence 
and ballistic reflection from the deeper lying atoms cause 
surface atoms to pop out. 
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FIG. 5. Projected trajectories for the adatorn/ vacancy 
forming process in the second column of Figure ^ (one of 
the 35 Ag — > Ag(lll) impacts). The horizontal coordinate is 
the one along the atomic row of atoms 3, 2, 4 and 5 in their 
initial position. We see atom 1 coming in and hitting 2 and 
4 (dashed lines), knocking them down and to the sides. 2 
and 4 bounce off the second layer and, still moving outwards, 
they push 3 and 5 out of the surface and take their places. 
The dots show the atomic positions after 0.25 ps. This is an 
example of a 35 eV trajectory, but it is also by far the most 
frequent mechanism of forming adatom vacancy pairs at 25 
eV. How far atoms 2 and 4 move down may vary — the main 
feature is that the incoming atom pushes two surface atoms 
apart, and these then each push a neighbor out of the surface 
layer. 
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FIG. 6. Average number of adatom-vacancy pairs created 
in Ag— >Ag(lll) impacts on an initially flat surface at 200 K as 
a function of the impact energy. Each data point is based on 
200-300 hundred simulated impacts, and the error bars give 
the statistical uncertainty. The 3 setups S2-S4 are deflned in 
Figure bI 



Figure || shows the average number of adatom/ vacancy 
pairs formed due to energetic Ag— >Ag(lll) impacts, as 
a function of energy, and for the three system setups 
S2, S3 and S4 discussed above. The production of 
adatom/ vacancy pairs is too high in the relatively small 
system with no Langevin atoms 5*2, and we have dis- 
cussed above how this is due to energy reflected from 
the system boundaries coming back to the local impact 
region. However, by introducing the shell of Langevin 
atoms as in system 6*3, we find an adatom/ vacancy pro- 
duction which agrees within the error bars with that of 
the significantly bigger system S'4. Based on this, unless 
otherwise stated, in the following we use systems at least 
as large as S3, i.e., iV = (> 4, > 5,3,4), with Langevin 
coefficients ^ = (0.005,0.010,0.15). Ni and N2 can be 
increased at higher energies. 

Figure ^ shows a significant produc- 
tion of adatom/ vacancy pairs for Ag— >Ag(lll) impacts 
with incoming energies greater than 20 eV. We shall see 
how this will influence the growth of the surface in sec- 
tion 



III B| . We have tested how this adatom/ vacancy pro- 



duction changes when we change the angle of incidence 
of the energetic beam. For 20 eV and 25 eV Ag/Ag(lll) 
impacts, at an angle 30 degree off normal incidence, we 
find average adatom-vacancy productions per impact of 
0.1 ±0.02 and 0.23 ±0.03, respectively. This calculation 
was done in a large cell with Ni = 8,7V2 = 6. Thus, at 
this angle, the adatom/ vacancy production is essentially 
the same as that for normal incidence. 
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FIG. 7. Possible outcomes of 25 eV impacts above a 
straight step on Ag(lll) initially at 200 K. From the left the 
columns are top view, initial and final configuration, then a 
view from an angle, initial and final configuration. The atoms 
which shift lattice sites are numbered, with 1 as the incom- 
ing atom. (The distance to the step is between -1 and -2 on 
Fig. ^.) From the top row down we see examples of 1) the 
equivalent of adatom/vacancy formation resulting in a dimer 
above the step; 2) no interlayer mobility, but a change of the 
step structure; 3) and 4) a net downward mobility by inser- 
tion of the incoming atom, as well as a change of the step 
structure. 2) and 4) also show the breaking off of adatoms 
from the step. 



resulting adatoms on the upper terrace have met to form 
a dimer in the first few picoseconds after the impact. 
The second row shows an example of no net interlayer 
mobihty due to the impact. The incoming atom gets in- 
corporated into the upper terrace, but another atom (la- 
beled 7) pops up. However, the impact does change the 
step structure, which is no longer straight. The collision 
has produced kinks on the step, and the atom labeled 3 
ends up detached from the step. This can influence the 
growth, because it might change the subsequent thermal 
diffusion. For example, the thermal interlayer mobility 
might be different at steps with kinks ||2^. 

Figure |^ rows 3 and 4 show examples of an energy- 
induced net downward interlayer mobility. Where for 
thermal impacts the incoming atom would stay on top 
of the upper terrace, it now gets inserted into the upper 
terrace due to the energy. This kind of process has pre- 
viously been seen in molecular dynamics simulations of 
Pt-^Pt(lll) impacts at similar energies [|l6|. Imagine the 
step is surrounding an island on the surface. The effect of 
the energy in this case is to change vertical growth of the 
island to horizontal growth at its periphery, and hence 
this insertion process will favor smooth growth of the 
surface. The examples show that the insertion process 
also causes a change of the step structure, and possibly 
detachment of atoms from the step to the lower terrace. 



3. Atom insertion 



We now proceed to examine the kind of energy- 
induced atomic rearrangements found for impacts near 
straight steps on the surface. We find it important to 
look at possible outcomes of these single impacts in or- 
der to identify and categorize the energy-induced atomic 
processes influencing the growth. On the other hand, it 
will soon become evident that it is not feasible to quan- 
titatively describe all possible outcomes of impacts on 
all possible atomic configurations of the surface. To see 
how energetic impacts may influence the total interlayer 
mobility during the growth we look at impacts in the 
vicinity of one of the two possible straight steps on the 
fcc(lll) surface: the so-called B-step exposing a (111) 
micro facet at the step edge. 

Figure ^ shows a selection of possible outcomes of 25 
eV Ag^Ag(lll) impacts just above a straight B-step. 
The flrst row shows the production of an adatom/vacancy 
pair, equivalent to what can take place on the flat sur- 
face, as discussed above. In this particular case the two 
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FIG. 8. Possible outcomes of 25 eV impacts below a 
straight step on Ag(lll) initially at 200 K. From the left 
the columns are top view, initial and final configuration, then 
a view from an angle, initial and final configuration. The 
atoms shifting lattice site are numbered, with 1 as the in- 
coming atom. (The distance to the step is between 1 and 2 
on Fig. ^.) From the top row down we see examples of: 1) 
an exchange event with no net change of the substrate; 2) 
an adatom/ vacancy formation resulting in attachment to the 
step; and 3) and 4) net upward mobility where a step atom 
gets piled up on the island. 



o 

Q. 

E 

1?) 
E 
o 



J2 
o 
E 

T3 
CD 
O 

T3 



1.0 
0.5 
0.0 
0.5 
0.0 
0.5 
0.0 
0.5 
0.0 




25 eV 



a h S' < A \ A ' A-" \ 



18 eV 



11 eV 



-4-2 2 
Distance to step 



4- Atom 



-up 



FIG. 9. Impact-induced interlayer mobility as a func- 
tion of the distance to a straight B-step for Ag^Ag(lll) 
(200K) at different energies. We show the average number of 
adatom/ vacancies formed (filled circles. Fig. ^ row 1, Fig. ^ 
row 2); the probability of insertion events (open squares, 
Fig. ^ rows 3 and 4) ; and the average number of atoms piled 
up on the island (open triangles. Fig. ^ rows 3 and 4). The 
"distance to step" is the horizontal distance in lattice con- 
stants to to the center of the step edge atoms (as shown at 
the top). We average over impact parameters within the same 
atomic row parallel to the step, and plot the result at the half 
integer distance (impact parameters between and 1 are av- 
eraged over, and plotted at 0.5 etc.). Statistical errors are less 
than 4% . 



Figure ^ shows examples of possible outcomes of 25 eV 
Ag^Ag(lll) impacts just below a straight B-step. The 
first row shows an impact which results in an exchange 
process involving two surface atoms and no net interlayer 
mobility. The only net effect of the energy is a short 
ranged horizontal displacement of the adatom from the 
impact site, as we have seen also happens for impacts 
on flat parts of the surface. The second row shows an 
example of an adatom/ vacancy formation on the lower 
terrace. This process is the same as the one for which 
the atomic trajectories are shown in figure |5[ only here 
the two adatoms can attach to the ascending step. 

Figure ^ rows 3 and 4 show examples of an energy 
induced net upward interlayer mobility. The atomic tra- 
jectories for the process in row 3 are very similar to those 
for row 1, only now atom 3 pops up in the step edge, and 
the step atom 4 gets piled up on the upper terrace. If 
the step is surrounding an island, this pile-up process 
causes vertical growth of the island, and hence it favors 
a rougher growth. 



5. Dependence on energy and step position 

Let us now turn to a quantitative discussion of the 
probabilities of the various interlayer processes we have 
seen can happen for impacts near straight steps. Again 
we look at impacts near the straight B-step, and em- 
phasize that this is merely one example of many surface 
configurations an incoming energetic atom can encounter. 
Figure |^ shows the average net number of atoms per im- 
pact moved up or down by adatom/ vacancy formation, 
insertion, or pile-up as a function of the distance to the 
step, and for 4 different incoming energies in the case 
of Ag— >Ag(lll). At 11 eV the only significant energy- 
induced process is the insertion mechanism. For impacts 
between the first and the second row of atoms in the up- 
per terrace, there is roughly a 50% chance of having an 
insertion event, hence there is a downward mobility of 
roughly 0.5 atom per impact. For impacts between the 
second and the third row of atoms, this probability has 
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decreased to approximately 10% . There is a few percent 
chance of insertion events for impacts between the third 
and the fourth atomic row in the upper terrace. Pile-ups 
are defined to be when an adatom ends up on the upper 
terrace for impacts outside the crystal atomic position of 
the step atoms, i.e., for "distance to step" greater than 
zero. At 11 eV we see that there is a few percent chance 
of pile-ups for impacts within one atomic row width from 
the step. We note that with this definition, it is possible 
that there would be a few pile-ups even thermally for 
impacts very close to the step edge. 

At 18 eV, the picture is more or less as at 11 eV, except 
that there is an increased probability for the insertion 
events, in particular for impacts between the second and 
fourth atomic row in the upper terrace. We also see a 
few percent chance of adatom/ vacancy production. The 
increase in the insertion events does not continue with 
increasing incoming energy. At 25 eV it is slightly re- 
duced compared 18 eV, and at 32 eV is further reduced 
as adatom/ vacancy production becomes more dominant. 

But the number of pile-up events has increased with 
incoming energy, and at 32 eV, there are more pile-up 
events for impacts below the straight step than inser- 
tions for impacts above the step. Both pile-ups and in- 
sertions happen within the width of three to four atomic 
rows from the step. In agreement with figure |^ the av- 
erage number of adatom/ vacancy pairs per impact has 
increased at 25 and 32 eV to approximately 0.3 and 0.5 
for impacts not too close to the step. However, there is a 
clear dip in the probability for these events close to the 
step, in particular for impacts right below the step. 

Based on figure ^, one might conclude that growing a 
Ag(lll) surface using a low energy atom or ion-beam 
might have an advantageous effect on the surface rough- 
ness, since at 18 eV we see a significant energy-induced 
downward mobility near step edges. However, one would 
also expect an optimal energy for growing a smooth sur- 
face, since at a little higher energy the total pile-up ex- 
ceeds the total insertion at the straig ht step. This is in- 
deed what we will see in section EH B in a certain growth 
parameter range. However, it must be kept in mind that 
the importance of insertions and pile-ups will depend 
on the step density on the surface, and that when the 
steps on the surface are rough the probabilities for inser- 
tions and pile-ups might be different from those for the 
straight step. 



at these energies. At 25 eV there is a few percent chance 
of producing adatom/ vacancy pairs during the impact 
and, by comparison to Ag— >Ag(lll), we conclude that 
this process must have an onset at a higher energy. 
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FIG. 10. Impact-induced interlayer mobility as a func- 
tion of the distance to a straight B-step for Pt— >Pt(lll) 
(200K) at different energies. We show the average number 
of adatom/ vacancies formed (filled circles); the probability of 
insertion events (open squares); and the average number of 
atoms piled up on the island (open triangles) — see text to 
Figure bl 




6. Pt-^Pt(lll) impacts 

Figure |l^ shows the results of the same kind of calcu- 
lation presented in figure ^, but now in the case of 18 eV 
and 25 eV Pt^Pt(lll) impacts. The energy-induced 
processes in the vicinity of the step are predominantly 
insertions. Both at 18 eV and 25 eV, there is a signif- 
icant probability of insertion events up to three to four 
atomic rows from the step. There are no pile-up events 



FIG. 11. Results of 18 eV Ag^Ag(lll) impacts on a small 
island initially with an adatom atop. The initial configuration 
is shown in the upper left — the incoming atom (labeled 1) 
is a few atomic distances above the island. Three possible 
outcomes are: 1) two atoms atop the island (upper right); 2) 
one atom atop the island (lower left); and 3) no atoms atop 
the island (lower right). 
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7. Insertion of additional atoms and island break-up 

Another mechanism favoring smooth growth was re- 
vealed to us watching movies of the full growth from the 
KMC-MD simulations. In impacts on a small island with 
an adatom atop, it is possible that the incoming atom and 
the adatom initially atop the island both end up in the 
layer of growing island. Thus, not only is the incoming 
atom inserted into the island, an existing adatom on the 
island is inserted as well. In these processes the small 
island is likely to undergo structural changes, perhaps 
even break up. We tried thirty 18 eV Ag^Ag(lll) im- 
pacts on a 3 by 3 atom island with an adatom initially 
on top, with random impact parameters in a surface unit 
cell neighboring the top adatom (see fig. Of these 

30 impacts, both atoms ended up in the growing layer 
in 19 cases (2 atoms down) , one atom stayed on top of 
the island in 7 cases (1 atom down), and in 3 cases both 
atoms stayed on the island (no induced interlayer mobil- 
ity). In one case 3 atoms ended up on top of the island 
(1 atom up). Regarding the impact induced break-up 
of the island, only in 8 cases did the 11 atoms end up 
as one connected island, while in 22 cases one or more 
atoms ended up separated from the original island. 

8. Summary of MD simulations of single impacts 

In summary we can categorize the collision-induced 
events as follows: 

• Adatom/ vacancy formations. The impact 
causes the digging up of atoms from the surface 
layer. 

• Insertions. Near descending steps, the incoming 
atom can be incorporated into the upper terrace. 
Insertion can also happen for impacts on small is- 
lands. Adatoms initially residing above descending 
steps can also be inserted due to an energetic im- 
pact. 

• Pile— ups. Near ascending steps, impacts below 
the step can result in net growth on the upper ter- 
race. 

• Step edge restructurings. Impacts near steps 
can change the step edge structure. 

• Island break-ups. Impacts on small islands can 
cause the island to break up into several pieces. 
Also, impacts near straight steps can result in 
adatoms detaching from the step and ending up 
isolated. 

• Horizontal mobility. Incoming atoms may be 
displaced horizontally before thermalizing. Also, 
all the events above involve some induced horizon- 
tal mobility. 



Some of these events must be characterized as energy- 
induced defect formations, leaving the surface in a higher 
energy state when compared to a thermal deposition. 
Adatom/ vacancy formations, pile-ups and island break- 
ups are such examples. Insertions, on the other hand, can 
take the surface to a lower energy state, by incorporat- 
ing the incoming atom or other adatoms initially resid- 
ing above descending steps. In this respect insertions are 
examples of energy-induced annealing of defects. How- 
ever, as we have seen, insertions may be accompanied 
by an energy increasing restructuring of the step edge, 
by creating step edge defects like kinks. Also, step edge 
restructurings are defect formations for impacts on ini- 
tially straight steps, but for impacts near rough steps, 
they could be annealing events. 

One of the questions one might hope to answer by do- 
ing these molecular dynamics simulations of single im- 
pacts is what the nature of the energy-induced mobility 
is. Does the energy, for example, cause a local heating, 
and can we then understand the increased mobility to 
be a result of a locally higher temperature? Or is the 
impact induced mobility more ballistic in nature? Based 
on the results presented here, we conclude that for these 
low-energy (< 32eV) metal-on-metal impacts, the local 
heating and increased temperature is the wrong picture. 
Plots like the one shown in figure H] indicate that the more 
relevant picture is one of a clearly ballistic collision se- 
quence that, depending on impact parameter, may result 
in defect formation or annihilation or both, as discussed 
above. 



B. KMC-MD simulations of growth 

After having discussed what can happen in the ener- 
getic atom-surface collisions by presenting our molecular 
dynamics simulations of such impacts, let us now turn to 
the results of the combined KMC-MD simulation of the 
growth by energetic deposition. We begin by looking at 
the effect of the energy on the submonolayer growth of 
Ag/Ag(lll). Subsequently we discuss the growth of sev- 
eral layers, and then we discuss how the picture changes 
for Pt/Pt(lll). All the presented KMC-MD simulations 
are done at a deposition rate of 1 ML/s and with normal 
incidence. 
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FIG. 12. The surface morphology at 0.3 ML coverage for Ag/Ag(lll) as obtained for thermal growth and for growth by 
energetic deposition at four different incoming energies. The surface areas with periodic boundary conditions, 30x30, 40x40, 
and 50x50 atomic surface unit cells (decreasing with increasing energy), are chosen to give a minimum of 4-5 islands in the 
simulation. First and second layer growth is seen in all cases. The surface temperature is 60 K. 
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1. Ag/Ag(lll) submonolayer structure 

We have simulated Ag/Ag(lll) growth at four differ- 
ent incoming energies, 11, 18, 25 and 32 eV, and for 
thermal deposition, in which case no MD simulations 
are done. For all the simulations of energetic deposi- 
tion, we simulate the atom-surface collisions using molec- 
ular dyn amics for cluster setups as explained in sec- 



tion II C . All cluster setups have Langevin coefficients 
C = (0.005,0.010,0.15). At 11 eV we use TV = (3,4,3,4), 
at 18 eV and 25 eV we use N = (4,5,3,4), and at 32 
eV we use N — (5,6,3,4). We first show results where 
the surface temperature is set to 60 K. This temperature 
determines the rates of various diffusion processes in the 
KMC part of the simulation, and is also the Langevin 
temperature and the temperature to which the MD- 
atoms are initially equilibrated. We later discuss the 
effect of increasing the surface temperature. 

Figure [T2| shows the surface morphology at 0.3 ML cov- 
erage for Ag/Ag(lll) at the different energies, and for 
the thermal deposition. Focusing on the thermal run, we 
observe island formation due to surface diffusion. These 
islands are dendritic with very irregular step edges due to 
limited diffusion along the step edges at this low temper- 
ature. We also observe that the branches of the dendritic 
islands have preferred growth directions perpendicular to 
the so-called A-steps. This is because of the asymmetry 
in the corner diffusion in the EMT barriers, see table |. 
When a diffusing adatom attaches to the island at the 
abundant corner-sites, they subsequently preferenti ally 
move to the A-step, causing growth in that direction P6| . 

When comparing the thermal run with the various en- 
ergetic depositions in figure there is one striking dif- 
ference. For the energetic depositions, in particular at 25 
and 32 eV, we observe a higher density of smaller islands. 
This is even more clearly shown in figure |l^, where we 
show the island densities as a function of coverage in the 
lower panel. The figure shows that for not too high cov- 
erages the island densities during energetic depositions 
are higher than that for the thermal deposition — for 25 
and 32 eV they are significantly higher. The mechanism 
giving rise to this increased island density at the higher 
energies is the formation of adatom/ vacancy pairs in 
the energe tic atom-surface collisions, as discussed in sec- 



tion 
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Dimers are either created directly, or they are 
formed with an increased probability by the two adatoms 
that result from a collision-induced adatom/ vacancy cre- 
ation. At 32 eV even trimers may be formed directly. 
The formation of a dimer does not necessarily nucleate a 
new island, since these dimers are mobile in the KMC- 
diffusion model, and hence can diffuse on the surface and 
attach to an existing island. However, the dimers move 
more slowly than single atoms — they have a shorter 
diffusion length on the time-scale of the deposition rate, 
resulting in a higher island density. 
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FIG. 13. Island densities (lower panel) and step densities 
(upper panel), for Ag/Ag(lll) at various impact energies and 
thermal deposition for submonolayer coverages. The surface 
temperature is 60 K. At 25 and 32 eV the island density is 
significantly increased compared to thermal deposition due 
to adatom-vacancy formation in impacts. This results in an 
increased step density for these energies, enhancing the to- 
tal energy-induced insertion at steps. At 11 and 18 eV the 
island density is increased slightly above the thermal island 
density due to the breaking up of existing islands in impacts. 
Islands formed this way coalesce quickly, and the mechanism 
does not increase the step density. After 0.5 layers down, the 
step densities for the energetic depositions level off due to 
the smooth growth, whereas the thermal step density contin- 
uously increases as the surface roughens. 



In section III A we found the energy onset of collision 
induced adatom/ vacancy formations to be around 20 eV. 
The lower panel of figure ^ shows an increased island 
density even at 11 eV. At this energy the island den- 
sity is increased around 0.1 ML due to energy-induced 
breaking up of existing islands into smaller pieces. Fig- 
ure ^4^-d shows the evolution of an island during the 
11 eV energetic deposition. From a to b, the island has 
grown from 17 atoms to a total of 22 atoms, but is also 
broken into two pieces. The breaking up of the island is 
not due to thermally activated mobility, but rather is the 
result of an impact on the very narrow island. From b to 
c the island grows further, and then from c to d one of 
the pieces again breaks into two. This break up is also 
caused by an energetic impact. The island now consists 
of three pieces. And it is counted as three when the is- 



land density plotted in figure 13 is calculated. From d to 
e two of the pieces have coalesced by further growth, and 
in f at a total coverage of 0.15 ML, the island has grown 
and coalesced into one connected piece. 

For thermal deposition the island density is a very use- 
ful concept, because the average distance between islands 
actually is a good measure for the typical island separa- 
tion. The reason for this is that the islands are nucleated 
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by the random diffusion and aggregation of adatoms. In 
the vicinity of existing islands, adatoms diffuse and at- 
tach at the island edge, and this gives rise to a "denuded 
zone" around the island with a low adatom concentra- 
tion, and hence with strongly reduced probability of nu- 
cleating new islands. As a result islands on the surface 
tend to be approximately equally spaced, with a spac- 
ing determined by the diffusion of the adatoms. As also 
discussed by Esch et al. Q, this does not necessarily 
hold for energetic deposition, since islands can be cre- 
ated very close to existing ones. Here it happens by the 
breaking off of pieces of existing islands. It could also 
possibly happen by the direct nucleation of new islands 
by impacts on the surface near existing islands. 

In the thermal case, because of the separation of the is- 
lands, coalescence sets in at 0.4-0.5 ML coverage, and the 
island density then drops down. The example in figure |lj 
of energetic deposition shows coalescence of the pieces of 
the island as early as 0.1 ML. The little islands created by 
breaking off pieces of existing islands do not survive for 
long. In the upper part of the pictures in figure [ija-f is a 
separate island which has been nucleated independently 
by the diffusion of adatoms. The two islands in figure |ljf 
are both nucleated by surface diffusion and behave much 
like islands formed in thermal deposition. They coalesce 
around 0.5 ML coverage. 




FIG. 14. Evolution of an island at low coverages for 11 eV 
Ag/Ag(lll) deposition at 60 K. a-f are at coverages of 0.056, 
0.063, 0.081, 0.088, 0.131, and 0.150 ML respectively. We see 
that during the deposition and the growth of the island it first 
breaks into two and then into three pieces, which coalesce 
again before a coverage of 0.15 ML. Such early coalescence 
can be expected for islands formed by the beam breaking off 
pieces of an existing island, and is generally not expected for 
thermal deposition. At the top of the images, another island 
is growing. The islands nucleated by diffusion on the terrace 
are well separated, and will typically coalesce around 0.5 ML 
coverage or later. 



bigger island. At 11 and 18 eV the island density plotted 
in the lower panel of figure ^ can be interpreted as the 
thermal island density plus short-lived fluctuations above 
this level due to breaking off pieces of islands. 

At 25 and 32 eV adatom/ vacancy pairs are formed in 
the atom-surface collisions. Here the island density is 
determined first by the fraction of impacts resulting in 1 
adatom, 2 adatoms, dimers etc. and then the diffusion of 
these species. In addition to this, there will be relatively 
short lived fluctuations due to the breaking up of existing 
islands into pieces . 

In section 



[II A we identified energy-induced processes 



taking place near steps on the surface which could affect 
the smoothness of the grown surface because they involve 
interlayer atomic mobility. The importance of these pro- 
cesses is of course determined by the total step density 
on the surface during the growth. In the upper panel of 
figure ^ we plot the step density, defined as the number 
of lattice sites occupied by atoms with less than 6 neigh- 
bors in the same layer divided by the number of sites in 
one layer (the unit is ML — monolayers). For thermal 
deposition the step density can be estimated from the is- 
land density, the coverage and a fractal dimension of the 
island. For the energetic depositions, we see something 
different. The increase in island density below 0.5 ML 
at 11 and 18 eV, is not reflected in an increase in step 
density. The islands formed by breaking off pieces from 
existing islands do not contribute to the step density in 
the same way as islands nucleated by diffusion. 

At 25 and 32 eV, on the other hand, we see a sig- 
nificantly increased step density at not too high cover- 
ages. This is due to islands nucleated from the relatively 
slow diffusion of dimers resulting from collision-induced 
adatom/ vacancy creations. But again, we do not have a 
simple relation between island density and step density. 
At 25 eV we observe higher island density than at 32 eV, 
but we observe higher step density at 32 eV. While the 
step density is obviously an important quantity in ener- 
getic deposition, we conclude that the island density here 
is a much less useful concept than it has proven to be for 
thermal deposition. 

In summary, we find that energetic deposition in- 
creases island densities and step densities. Island densi- 
ties are increased by breaking existing islands into pieces, 
and at the higher energies because of collision-induced 
adatom/ vacancy formation. Only islands formed by the 
latter mechanism give rise to an increased step density. 



Also at 18 eV we see an increased island density setting 
in around 0.25 ML, and by watching growth movies, we 
identify the mechanism to be breaking off little pieces of 
the rather dendritic islands. These little islands usually 
have short lifetimes before they again coalesce with the 
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FIG. 15. Simulated "anti-phase" intensity for Ag/Ag(lll) 
at 60 K for thermal deposition, and for energetic deposition at 
11, 18, 25 and 32 eV. Energetic deposition induces oscillations 
indicating (transient) layer-by-layer growth. The strongest 
and slowest decaying oscillations are observed at 25 eV. 



2. Ag/Ag(lll) — growth of the first few layers 

We now move on to a discussion of our results for the 
growth of several atomic layers of Ag on Ag(lll). Our 
primary focus will be on the smoothness of the growth, 
and how the smoothness is affected by the incoming en- 
ergy 



Figure 15 shows the normalized anti-phase intensity 
which would be measured in anti-phase scattering of He 
H, RHEED, [|6), LEED |^ or X-rays |§, for our ther- 
mal and energetic deposition simulations. The intensity 
I is calculated from the simulations as a function of time, 
using the expression / — (X]i^o(~-'^)*(^*~^(i+i)))^' where 
6i is the fractional coverage in the i'th layer. 



25 eV 




FIG. 16. Surface morphologies for Ag/Ag(lll) at 60 K af- 
ter a total deposition of 3 atomic layers, for thermal (size 
50x50) growth and energetic deposition at 18 eV (size 40x40) 
and 25 eV (size 30 x 30). 



At 60 K there is no thermally activated interlayer mo- 
bility. Hence, for the thermal deposition, adatoms land- 
ing on top of existing islands will stay there. Islands are 
nucleated on top of islands, and multi-layer growth re- 
sults. The anti-phase intensity / drops rapidly to zero. 
Figure |l5| shows that, when the surface is grown by en- 
ergetic deposition, oscillations in / can be induced. The 
damped oscillations in / correspond to a decaying layer- 
by-layer growth where the surface partly recovers its ini- 
tial flatness after each atomic layer is deposited. 
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FIG. 17. Plot of the coverage in each growing layer (layer # 
labeled at each graph) above the initial surface as a function 
of the total deposition for Ag/Ag(lll) at 60 K. Upper panel 
is for 25 eV, and lower is for thermal deposition. 



That the surface grown by energetic deposition is more 
smooth than that grown by thermal deposition can be 
seen more directly from figure |l6|, which shows the surface 
morphology after 3 monolayers deposited, with the gray 
scale indicating the height above the initial surface, for 
18 and 25 eV and for thermal deposition. The same effect 
is evident in figure which shows the coverage in each 
atomic layer as a function of the total amount deposited 
for 25 eV and for the thermal deposition. After thermal 
deposition of three monolayers, growth of the 6th layer 
sets in while the second layer is still not completed. At 
the same total coverage in the 25 eV deposition, the third 
layer is almost complete, there is some growth in the 
fourth layer, and only a fraction of a percent of the fifth 
layer is occupied by atoms. 

The fact that using energetic deposition promotes the 
growth of flat surfaces is consistent with our results 
for the energy-induced atomic mobilit ies in the atom- 

At 60 K 
is the 



surface collisions presented in section pi A 
thermal interlayer mobility can be neglected. It 
energy-induced insertions for impacts near descending 
steps which lead to the growth of more smooth surfaces 
— impact parameters where the incoming atom would 
stay on top of an existing island in thermal deposition 
but, due to the energy, now result in lateral growth of 
the island. 

Figure nS shows that the magnitude of the oscillations 
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of / gradually increase when the energy is increased from 
eV to 25 eV. But by increasing the energy further to 32 
eV the oscillations are again reduced. There is an optimal 
energy around 25 eV for growing smooth surfaces. This is 
exactly what is predicted from figure which shows the 
statistics of possible outcomes of impacts near straight 
steps. Figure ^ shows a maximum of insertion events at 
18 eV, pile-ups setting in at 25 eV and dominating over 
the insertions at 32 eV. Based on figure ^, one might have 
expected an optimum for growing smooth surfaces at 18 
eV, and one might ask why we find this optimum to be 
closer to 25 eV. Several things play a role: 1) figure ^ 
shows the result for straight B-steps, but the islands at 
this temperature are quite small, and have irregular step 
edges. This might shift the balance of the insertion events 
to pile-ups to have the optimum at higher energies, and 
2) the overall importance of the energy-induced events 
near the steps will increase with the step density and, as 
we have discussed above, the step density increases with 
increasing deposition energy after the onset of energy- 
induced adatom/ vacancy formations. At 25 eV, where 
the insertion events still dominate over the pile-ups at 
the straight steps, the total energy-induced downward 
mobility is increased due to the energy-induced increase 
in step density. 
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mobility along the island edges has increased, making is- 
lands more compact, and reducing the step density even 
further. Thus the total energy-induced downward mobil- 
ity near the step edges is reduced, and the surface grows 
rough more rapidly. It should be evident that the impor- 
tant parameter here is not the surface temperature itself, 
but rather the step density, and possibly the step edge 
structure. Due to the very low barrier for surface diffu- 
sion of adatoms on Ag(lll) (which is underestimated in 
the EMT), we have to go to the low temperatures used 
in our simulations to see an effect of the energy on the 
smoothness of the surface grown. 

As the temperature increases even further, the island 
and step densities decrease to such an extent that the 
effect of the incoming energy on the surface height distri- 
bution is negligible. However, there can still be a strong 
effect of the energy on the island density. For example, 
in (shghtly flawed @) KMC-MD simulations of 25 eV 
Ag(lll) growth at 200 K the energetic beam no longer 
changed the surface height distribution, even though the 
island density changed by a factor of up to 20 compared 
to thermal growth. However, we have found this to be 
strongly dependent on details of the implementation of 
the KMC diffusion model. First of all, the effect de- 
pends strongly on the mobility of the small clusters of 
atoms nucleated in energetic impacts. In the KMC model 
presented in this paper, these clusters are rather mobile 
due to periphery diffusion of their atoms, hence reduc- 
ing the energy effect on island densities. Another subtle 
but important factor is the efhciency with which diffus- 
ing dimers and trimers can fill in the isolated vacancies 
created in energetic impacts. In EMT this readily hap- 
pens. However in KMC models, where we do not allow 
this recombination, the concentration of vacancies may 
build up to an extent where it hinders the diffusion of 
the small clusters of atoms, and hence results in a higher 
island density (a point also made by Esch et al. 



FIG. 18. Simulated "anti-phase" intensity for Ag/Ag(lll) 
at two surface temperatures, 60 and 105 K. At the higher 
temperature the oscillations decay more rapidly, essentially 
due to a lower island density, a lower step density, and hence 
a lower total energy-induced downward mobility. 



3. Higher surface temperatures 

This takes us to the discussion of the effect of surface 
temperature during growth by energetic deposition. Fig- 
ure [18| shows that, for growth at 25 eV, if we increase the 
surface temperature from 60 K to the 105 K the magni- 
tude of the first oscillation in / is reduced approximately 
by a factor of 3. At 105 K the thermal interlayer mobility 
is still very low, but the diffusion of the adatoms on the 
terrace has sped up significantly, reducing the island den- 
sity, and hence the step density. In addition to this the 



^. Pt/Pt(lll) suhmonolayer structure 

We now present the results of KMC simulations of ther- 
mal and KMC-MD simulations of 18 eV and 25 eV ener- 
getic deposition of Ft/Ft (111) and compare these results 
to the case of Ag/Ag(lll). All the Ft KMC-MD simu- 
lation are done with a surface temperature of 80 K and 
a deposition rate of 1 ML/s. The MD part of the simu- 
lations is done using setups with Langevin coefficients 
^ = (0.005,0.010,0.15) and N ^ {Ni, N2, N3, N4) = 
(4,5,3,4). 
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FIG. 19. The surface morphology at 0.3 ML coverage for 
Pt/Pt(lll) as obtained for thermal growth and for growth 
by energetic depositions at 18 eV and 25 eV. The surface 
areas with periodic boundary conditions are 30x30 and 50x50 
atomic surface unit cells for energetic and thermal deposition, 
respectively. First and second layer growth is seen in all cases. 
The surface temperature is 80 K. 

Figure |l^ shows the surface morphology after the depo- 
sition of 0.3 ML of Pt on Pt(lll) for thermal deposition 
and for the two energies. Compared to the Ag/Ag(lll) 
simulations, there is a much higher thermal island den- 
sity. This is of course due to the higher energy barriers 
for surface diffusion on Pt. Because of the higher density 
of islands, they are smaller and tend to be less branched. 
Thus in the case of Pt, the submomolayer structure ob- 
tained by energetic depositions at 18 and 25 eV is very 
similar to the thermal submonolayer structure. This is 
also apparent from figure pO| , which shows the island and 
step densities, as figure ll^ did for Ag. At low cover- 
ages the island and step densities for the energetic de- 
positions follow the corresponding curves for the thermal 
deposition quite closely. If there is a difference at all, the 
energy slightly reduces the island density. This could be 
explained by short ranged energy-induced horizontal mo- 
bility. Around 0.7 ML the energetic step densities level 
off, indicating smooth growth, whereas the thermal one 
continues to increase as the surface roughens. We recall 
from figure |l^ that 25 eV is below the onset for energy- 
induced adatom/ vacancy formations for Pt^Pt(lll) im- 
pacts, consistent with the unchanged island density. Ap- 
parently, the breaking up of existing islands is also not 
pronounced here — probably as a consequence of the less 
branched island structure, and the stronger Pt-Pt bond 
energies. 
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FIG. 20. Island densities (lower panel) and step densities 
(upper panel), for Pt/Pt(lll) at 18 eV, 25 eV and thermal 
deposition for submonolayer coverages. The surface temper- 
ature is 80 K. In the case of Pt, 18 and 25 eV are both too 
low energy to change the island and the step densities at low 
coverages. 



5. Pt/Pt(lll) — growth of the first few layers 




Total deposition, [ML] 

FIG. 21. Simulated "anti-phase" intensity for Pt/Pt(lll) 
at 80 K for thermal deposition, and for energetic deposition 
at 18 and 25 eV. Energetic deposition induces oscillations, 
indicating (transient) layer-by-layer growth. Fairly strong 
and slowly decaying oscillations are observed at 18 and 25 
eV. 



While the 18 eV and 25 eV energetic deposition has 
no influence on the surface morphology at low coverages, 
it has a significant influence as the growth progresses. 
Figure |2l| shows fairly strong and slowly damped oscilla- 
tions in the simulated anti-phase intensity being induced 
by the energy. At 25 eV the periodicity in the growth of 
each layer is evident from figure ^ {i.e.. the curves for 



18 



successive layers look the same) . For thermal growth, on 
the other hand, the completion of each layer spreads out 
over an increasing period of time. 
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FIG. 22. Plot of the coverage in each growing layer (layer ^ 
labeled at each graph) above the initial surface as a function 
of the total deposition for Pt/Pt(lll). Upper panel is for 25 
eV, and lower panel is for thermal deposition. 



Two main factors play a role in making the energetic 
growth so smooth in the case of Ft. First, figure [lO| 
shows a pronounced energy-induced mobility for impacts 
near steps that favors smooth growth, i.e., for impacts 
above the step, there is a high probability of energy- 
induced insertions. We do not see induced pile-ups for 
impacts below the step at these energies. Second, the 
step density for Ft/Pt(lll) at this temperature is very 
high. As was the case for Ag we expect the advantageous 
effect of the energy to die away with increasing surface 
temperature, as the step density decreases with increased 
surface diffusion. 
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FIG. 23. Surface morphologies for Pt/Pt(lll) at 7 meV 
(80 K) after a total deposition of 5 atomic layers, for thermal 
(size 50 X 50) growth and energetic deposition at 18 (size 
30 X 30) and 25 eV (size 30 x 30). 



IV. DISCUSSION OF RESULTS 

Here we summarize and discuss our results from simu- 
lating growth by energetic deposition. Our results from 
our study of the possible outcomes of single impacts on 
diffe rent su rface configurations were summarized in sec- 
tion III A 8 . We found energy-induced defect formation 
mechanisms, such as adatom/ vacancy formations, pile- 
ups, step edge restructurings and break-ups of existing 
islands. We also found energy-induced defect annihi- 
lation mechanisms, in particular insertion of atoms into 
horizontally growing layers. For the energies in this study 
we in general found the energy-induced mobility to be 
limited to a range of a few (~5) atomic distances. 

In doing the KMC-MD simulations of the entire growth 
process, we have seen that the energy-induced atomic 
mobility can affect island densities in two ways. First, 
for energies above a certain threshold (20 eV for Ag) , the 
energy-induced formation of adatom/ vacancy pairs and 
direct formation of dimers and trimers on the surface in- 
creases the island density. In our simulation these small 
clusters of atoms diffuse only slowly on the surface, giving 
rise to the reduced island separation. It is an important 
characteristic for islands formed in this manner that they 
are nucleated by diffusion and aggregation: they tend to 
repel each other via their diffusional fields, and are well 
separated as is the case for islands formed in thermal 
growth. Hence, they coalesce in the late stages of the 
growth of a monolayer. Second, the energetic impacts 
increase the island density in a separate way, by break- 
ing off clusters of atoms from existing islands. While 
this process may give rise to new islands, they are always 
nucleated very close to existing islands, and are qualita- 
tively different. In the Ag/Ag(lll) simulations, we have 
seen that they do not contribute to the step density, in 
contrast to islands nucleated by diffusion. They have a 
very short life time, either because they diffuse to join 
the original island, or because they grow and coalesce at 
a much earlier growth stage. 

We have seen that at low temperatures using energetic 
deposition can change rough thermal growth to smooth 
layer-by-layer growth. As usual, this is expected eventu- 
ally to decay into rough growth as the number of layers 
increases |Q. We have done simulations at tempera- 
tures where thermal inter layer mobility is frozen out, and 
have found that energy-induced insertions of atoms into 
growing layers can be sufficient to give layer-by-layer 
growth. However, energy-induced pile-ups, which set in 
at a higher energy, can cause the opposite effect. For 
this reason we find that there is an optimal energy for 
layer-by-layer growth — for Ag it is approximately 25 
eV. We expect the existence of such an optimum to be 
general. Pile-ups are likely to dominate over insertions 
at higher energies; for insertions only the incoming atom 
or other atoms residing above descending steps can be 
inserted, while the number of atoms that can pile up is 
not similarly limited. 
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Since the energy-induced insertions are limited to im- 
pacts within a short range of descending steps, we have 
seen that a high step density is needed for energy-induced 
smooth growth. While the use of energetic deposition can 
itself assist in increasing the step density, we have seen 
the energy induced layer-by-layer growth dies away with 
increasing temperature in the case of Ag/Ag(lll). The 
horizontal diffusivity for this system is very high, and the 
two dimensional islands readily grow very large as the 
temperature is increased, giving a very low step density. 
However, how efficient the energetic beam is in increasing 
the island density depends on many details of the atomic 
potential energy landscape. In some cases this efficiency 
can be very high. Therefore, in some systems which have 
a low step density when thermally grown, it is possible 
that an energetic beam might increase the step density 
sufficiently for the energetic-beam insertions to give rise 
to smooth growth. 



V. DISCUSSION OF METHOD 

We now turn to a general discussion of our Kinetic 
Monte-Carlo Molecular Dynamics method for simulat- 
ing crystal growth by energetic deposition. We wish to 
address the validity of the method, its advantages and 
disadvantages, and how we find it useful in providing in- 
formation about the energetic growth. 

Kinetic Monte-Carlo treats the thermal diffusion as 
a sequence of uncorrelated atomic hops, and assumes 
these diffusion hops to be instantaneous. In reality, a 
diffusion hop in a crystalline environment has a dura- 
tion time Tjn of roughly a picosecond, the typical time 
an atom spends on top of the barrier while crossing it. 
The Kinetic Monte-Carlo formalism is correct when the 
times between various kinds of hops is much larger than 
Tm- (The time between hops is the inverse rate of the 
diffusion process, given by the Arrhenius formula with 
the barriers in Table ^) In the same way, KMC-MD is 
correct when the time between hops is larger than te, 
the duration of our MD simulation of the non-thermal 
collision-induced mobility. During te we neglect ther- 
mal mobility outside and across the boundaries of the 
region of the MD simulation. This approximation only 
limits the validity of the simulation to the extent that 
this mobility would correlate to the atomic rearrange- 
ments in the impact region. The collision is treated as 
instantaneous by the KMC algorithm, so the total mobil- 
ity outside the impact region is correctly accounted for, 
whereas inside the region thermal mobility within te is 
doubly accounted for. While this should be kept in mind, 
we don't expect it to have any significance for the sim- 
ulations presented in this paper, because diffusion is so 
slow compared to te- 

The KMC-MD method also carries over the other 
usual limitations of KMC. One of these is that the KMC 
part of the simulation must be done on a lattice. In the 



present simulations of growth of fcc(lll) surfaces, we ne- 
glect the presence of the off-lattice hep binding sites in 
the KMC part of the simulation, and hence we have to 
move atoms trapped in these sites at the end of an MD 
simulation to neighboring fee sites. This is severe only 
if growth in practice would take place on the hep sites. 
Growth partly on hep sites, partly on fee sites, would give 
rise to complicated structures and dynamics at the stage 
of coalescence of these islands, a scenario we are unable 
to study with the KMC-MD method. As noted earlier, 
the favorable hep site is an artifact of the EMT poten- 
tial we use, and these issues are most likely not relevant 
exp er iment ally. 

The KMC-MD method must be based on a model 
potential for the atomic interactions, in this case the 
EMT. While the EMT includes many qualitative and to 
some extent quantitative features of the interaction of 
late transition and noble metals, it is not an exact po- 
tential. In the case of Pt(lll), for example, it is found 
that a simple scaling of the EMT energy barriers by a 
factor of 1.6 as input for a KMC simulation of thermal 
growth gives a good agreement with experimental island 
densities, transition of fractal to compact islands, and 
the appearance of reentrant layer-by-layer growth at low 
temperatures p5| , ^j27| . A common scaling factor of all 
energy barriers would not be a severe discrepancy — it 
merely corresponds to a scaling of the temperature. 

While the approximations made in effective medium 
theory prevent us from controlled predictions of exactly 
how a real material would grow by energetic deposition, 
the method still includes the right ingredients to be a 
useful tool for our purposes: to examine how the use 
of energetic particles may influence the crystal growth; 
and to identify mechanisms and evaluate their relative 
importance, as has been done in the previous sections. 
We have demonstrated that the KMC-MD method can 
indeed be used to elucidate the important interplay be- 
tween energy-induced mobility and thermal surface dif- 
fusion. 

One could imagine an alternative way of doing these 
simulations of growth by energetic depositions. First 
make a complete table of all possible outcomes of single 
atom-surface collisions listed together with their relative 
statistical significance, and for all relevant local surface 
configurations. Then perform a KMC simulation, which 
for each impact would choose the resulting local config- 
uration from such a table. The present simulation study 
sho ws wh at an enormous task that would be. In sec- 
tion III A we considered impacts on the fiat surface and 



near one type of straight step, and found a dependence on 
the step position up to 5 atomic distances away. Imagine 
the table of collision events needed to correctl y acc ount 
for the evolution surfaces presented in section |III B ! 



A major advantage of the present method is that it 
allows one to evolve the surface by depositing energetic 
atoms without making any assumptions about the effect 
of the energy. From each simulation we can in princi- 
ple record exactly the contributions of various types of 



20 



energy-induced atomic mobilities. The method correctly 
convolves the distribution of possible outcomes of ener- 
getic impacts on a given local configuration with the dis- 
tribution of different local configurations during the crys- 
tal growth. The energy-induced microscopic mechanisms 
acting are not assumed, but on the contrary revealed by 
the method. One example of this is the insertion of atoms 
residing atop islands prior to impacts, as discussed in sec- 
tion 
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The main disadvantage of the method is that it is very 
time consuming computationally. The MD simulations 
are expensive, and a new simulation for every deposited 
atom is needed. Even for the relatively small surface ar- 
eas considered here, a thousand to a few thousands of 
MD simulations are needed per monolayer. In order for 
the method to be feasible, it is essential that the island 
density not be too low. Big islands require a larger sur- 
face area and hence a larger number of MD simulations 
per monolayer. Our simulations ran for weeks to months 
on reasonable 1997 workstations. 

Our method is not well suited for parallelization. The 
bulk of the CPU time is running the MD for the atomic 
impacts. Each MD simulation is small and parallelization 
by subdivision would give rise to a considerable commu- 
nication overhead. The current size of the simulation is 
set by the need to have an impact not be affected by the 
periodic boundary conditions: the impact and its further 
diffusive evolution must not feel its periodic images. Dis- 
tribution of the impacts among different processors would 
necessarily involve artificial simultaneous depositions in 
different areas of the surface. For these depositions to be 
unaffected by one another, it would seem that the same 
system-size per processor must be needed: parallelization 
does not allow more monolayers per month. 

Our new method, by combining MD for the impact 
with KMC for the diffusive motion between impacts, al- 
lows for realistic simulations of surfaces grown with en- 
ergetic beams at realistic growth rates of monolayers per 
second. Good simulations, however, take a week or more 
per monolayer deposited. 
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